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ABSTRACT 

The  geometric  and  kinematic  description  of  the  motion  of  an  elastic-plastic 
body  undergoing  finite  deformation  is  considered.  A  configuration  distinct  from 
the  initial  and  current  configurations  of  the  body  is  introduced,  which  is 
supposed  to  be  the  state  of  the  body  due  to  its  plastic  deformation  alone,  and 
also  to  be  the  reference  configuration  for  elastic  deformation.  It  is  shown 
that  the  rate  of  deformation  tensor  may  then  be  written  as  the  sum  of  the  rates 
of  elastic  and  plastic  deformation.  Restrictions  imposed  by  some  physical 
assumptions  are  considered,  and  the  form  of  some  of  the  derived  relations  is 
given  for  two-dimensional  deformations  in  Cartesian  and  cylindrical  coordinates. 
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-  deformation  tensors 

-  Green's  tensors 

-  stretching 
strain 

-  metric  tensor 

-  associated  metric  tensor 

-  Euclidean  shifters 

-  measure  of  dilation 

-  inverse  deformation  tensors 

-  polar  coordinates 

-  arc  length 

-  time 

-  material  coordinate 

-  reference  coordinate 

-  current  coordinate 

-  Cartesian  coordinates 

-  cylindrical  coordinate 

-  unit  tensors 

-  Christoffel  symbols 

-  relative  permutation  tensors  of  weight  +1  and  -1,  resp. 

-  cylindrical  coordinates 

-  cofactor  matrix  of  (  ) 

-  covariant  derivative 

-  determinant 

-  material  time  derivative 

-  physical  component  of  a  tensor 

-  symmetric  part  of  (  ) 

-  trace  of  (  ) 
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INTRODUCTION 

In  recent  years,  interest  in  problems  concerned  with  large  stress,  high 
speed  loading  of  bodies  of  elastic-plastic  material  has  been  generated  due  to 
several  important  applications.  Among  these  are  penetration  of  armor  plat¬ 
ing  by  projectiles  and  the  use  of  explosives  in  metal  forming.  Because  of  the 
large  traction  amplitudes  and  high  loading  rates  involved,  it  has  become 
necessary  to  construct  constitutive  models  for  such  materials  for  the  case  of 
finite  deformations.  Before  the  existing  theories  of  plasticity  can  be  gener¬ 
alized  to  the  case  of  finite  deformation,  however,  certain  questions  about  the 
geometric  and  kinematic  description  must  be  reconsidered.  For  example,  the 
most  obvious  restruction  imposed  by  introduction  of  finite  deformation  is  that 
the  elastic  and  plastic  strains  are  no  longer  simply  additive. 

Even  though  consideration  here  is  limited  to  discussion  of  kinematics  and 
geometry,  the  terminology  associated  with  a  particular  material,  namely 
elastic-plastic,  is  used.  This  is  because  the  work  here  is  intended  solely 
for  applications  involving  such  materials.  As  usual,  the  descriptions  elastic 
and  plastic  deformation  refer  to  the  recoverable  and  irrecoverable  parts, 
respectively,  of  a  homogeneous  deformation  of  an  elastic-plastic  body. 

As  is  well  known,  the  description  of  the  motion  of  an  elastic  or 
viscoelastic  body  may  be  given  by  considering  the  relation  between  two  config¬ 
urations  of  the  body,  the  initial  and  the  current  configurations.  To  discuss 
the  motion  of  an  elastic-plastic  body  it  is  found  to  be  convenient  to  introduce 
a  third,  intermediate  configuration.  Thus,  in  addition  to  the  initial  and  cur¬ 
rent  configurations, we  consider  a  varying  intermediate  configuration  which  is 
the  state  the  body  would  be  in  due  to  its  plastic  deformation  alone,  and  which 
is  treated  as  the  reference  configuration  for  elastic  deformation. 
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With  the  introduction  of  this  intermediate  state,  it  is  shown  that  the 
total  rate  of  deformation  or  stretching  tensor  can  be  written  as  the  sum  of 
the  rates  of  elastic  and  plastic  deformation.  The  significance  of  this  is 
recognized  when  one  considers  the  usual  development  of  constitutive  equations. 
An  appropriate  law  of  elasticity  gives  the  rate  of  elastic  deformation  in  terms 
of  stress  and  stress  rate.  In  addition,  the  rate  of  plastic  deformation  is 
prescribed  in  terms  of  stress  by  an  appropriate  flow  rule.  These  two  terms 
can  then  be  added  to  yield  the  total  rate  of  deformation  in  terms  of  stress. 

The  idea  of  introducing  an  intermediate  configuration  to  discuss 
elastic-plastic  deformation  was  apparently  first  put  forth  by  Eckart  [1]  to 
discuss  materials  with  rate  effects.  Related  discussions  concerned  with  kine¬ 
matics  or  particular  constitutive  laws  have  been  given  by  Backman  [2], 

Eglit  [3],  Sedov  [4],  Truesdell  [5],  and  Truesdell  and  Toupin  [6].  The  con¬ 
cepts  were  finally  put  to  use  by  Lee  and  his  colleagues,  as  discussed  in  a 
series  of  papers  on  the  propagation  of  large  amplitude,  one-dimensional  strain 
waves  [7,8,9]. 

The  object  here  is  to  generalize  the  kinematical  concepts  given  by  Lee 
et  al.  to  the  case  of  orthogonal  curvilinear  coordinates.  To  accomplish  this, 
use  is  made  of  the  standard  tensor  notation  and  the  theory  of  double  tensor 
fields,  as  presented  by  Ericksen  [10].  For  a  discussion  of  the  general  theory 
it  would,  of  course,  be  simpler  and  more  reasonable  to  consider  only  Cartesian 
coordinates.  We  have  particular  applications  in  mind,  however,  which  are  best 
suited  for  curvilinear  coordinate  systems.  The  slight  extra  effort  required 
here  will  hopefully  be  justified  when  the  results  are  applied  to  these 
particular  problems. 

After  the  general  geometric  and  kinematic  relations  have  been  derived,  the 
forms  of  the  assumptions  of  zero  plastic  volume  change  and  infinitesimal  shear 
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deformations  are  investigated.  Finally,  brief  remarks  are  made  with  regard  to 
two-dimensional  deformations  in  Cartesian  and  cylindrical  coordinate  systems. 


GEOMETRIC  AND  KINEMATIC  RELATIONS 

Consider  a  simply-connected  body  of  elastic-plastic  material,  prior  to 

application  of  any  loading,  whose  material  particles  may  be  located  by  a  fixed 

Cartesian  coordinate  system.  A  fixed  set  of  orthogonal  curvilinear  coordinates 
K 

Xq  is  then  prescribed  in  the  region  occupied  by  the  body,  and  each  material 
particle  is  identified  with  the  coordinate  Xq  which  it  occupies,  i.e.,  the 
Xq  are  material  coordinates.  Since,  in  this  configuration,  the  material  co¬ 
ordinate  system  possesses  an  underlying  Cartesian  coordinate  system  the  corre- 

1  K 

sponding  metric  space  is  Euclidean,  and  the  distance  dSQ  between  points  Xq 

and  X^  +  dX^  is  given  by 

dS2  =  G„r  dXK  dXL  (1) 

o  o  KL  o  o 

The  tensor  G„r  is  the  Euclidean  metric  tensor  of  the  curvilinear  coordinate 
o  KL 

system.  If  the  system  itself  is  rectilinear,  this  metric  is  just  the  identity 
tensor. 

Suppose  that  now  the  body  is  subjected  to  (time-dependent)  loading, 

resulting  in  a  change  of  configuration  and  temperature  distribution.  It  is 

assumed  that  the  loading  is  severe  enough  so  that  part,  and  perhaps  all,  of  the 

body  undergoes  some  plastic  deformation.  The  particle  which  was  originally  at 
K  i 

Xq  now  moves  to  the  place  x  .  The  deformation  carrying  this  arbitrary  parti- 
K  i 

cle  from  Xq  to  x  is  assumed  to  satisfy  the  axiom  of  continuity,  as  well  as 

■^The  word  "space"  is  used  here  merely  to  denote  the  totality  of  points 
corresponding  to  all  values  of  the  coordinates  within  certain  ranges.  This 
is  to  be  distinguished  from  "metric  space,"  which  is  a  space  in  which  the 
concept  of  length  has  been  introduced. 
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the  principles  of  impenetrability  and  permanence  of  matter.  The  deformation 
is  then  a  uniquely  invertible  topological  mapping,  and  is  described  by  the 
differentiable  function 

x1  =  x1(XK  ,  t)  ,  (2) 

o 

which  determines  the  current  configuration  of  the  body.  Explicit  dependence 
on  time  is  included  because  the  loading  may  be  time-dependent  or  the  material 
may  be  rate -dependent.  The  motion  (2)  prescribes  the  rule  not  only  for  the 

motion  of  an  arbitrary  particle ,  but  also  for  the  deformation  of  the  neighbor- 

K  K 

hood  of  the  particle.  Thus,  if  dXQ  is  a  line  element  emanating  from  Xq  , 

it  is  deformed  into  the  line  element  dx1  emanating  from  x1  by  the  differ¬ 
ential  relation 

dx1  =  x^K  dX*  =  a1KdX^  .  (3) 

The  tensor  a1^  is  the  deformation  gradient. 

Since  x1  is  also  a  point  in  the  aforementioned  space  of  orthogonal 
curvilinear  coordinates,  the  metric  tensor  g..,  which  determines  distance 

i] 

between  x1  and  x1  +  dx1  by 

ds2  =  g„  dx^  dx^  ,  (4) 

* 

has  the  same  components  as  except  that  they  are  evaluated  at  x 

K  K  i 

rather  than  at  Xq  ,  i.e.,  if  G^  =  F(Xo>,  then  g^  =  F(x  ).  The  mathemati¬ 
cal  relation  between  G.,r  and  g..  can  be  given  by  introduction  of  the 

o  KL  6i]  6  J 

K  i 

Euclidean  shifter.  If  Yq  and  y  are  the  points  in  the  underlying 

K  i 

Cartesian  system  corresponding  to  Xq  and  x  ,  then  there  exist  invertible 
coordinate  transformations  XK  =  XK(Y^)  and  x1  =  x1(y^).  The  shifters  g* 
and  g.^  (from  X^  to  x^")  and  g^.  and  gv^  (from  x^  to  X^)  are  then 

1  O  1  K  O 


defined  as 
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K 


1  =  si  Us.  /  . 

<  UT  n  V  »  &  J 


L  3y3  3X 


K 


3xK  i 
■\  _£  ill 

^  3YL  3xX 
o 


(5) 


•  ^ 

where  6  and  6  .  are  the  shifters  for  a  rectilinear  coordinate  system  and 
L  X 

equal  1  or  0  according  as  the  indices  are  the  same  or  different.  The  rela¬ 
tion  between  metrics  at  the  two  points  is  then 


_  K  L 
oGKL  g  ig  j  * 


(6) 


Relation  (6)  may  be  verified  by  direct  substitution. 

In  an  attempt  to  separate  the  elastic  deformation  from  the  plastic 
deformation,  we  follow  Lee  and  Liu  [7]  in  defining  a  particular  intermediate 
state.  This  configuration  is  the  one  that  the  body  would  be  in  due  to  its 
plastic  deformation  only,  that  is,  the  stresses  in  the  body  which  resulted  in 
elastic  deformation  are  completely  relieved  and  the  local  temperature  is 
reduced  to  its  initial  value.  For  a  body  which  has  undergone  nonhomogeneous 
plastic  deformation  this  configuration  cannot,  in  general,  be  achieved;  it  is 
a  conceptual  configuration  rather  than  a  physical  one.  This  is  due  to  the  fact 
that  when  the  loads  and  temperature  variations  are  removed  from  the  body,  after 
resulting  in  nonhomogeneous  plastic  deformation,  a  residual  stress  and  associ¬ 
ated  elastic  strain  remain  in  the  body.  To  reach  the  desired  stress-free 
configuration,  the  body  must  be  considered  to  be  cut  up  into  small  elements. 
When  a  material  element  is  removed  from  its  neighboring  elements,  which  have 
restrained  it,  the  stress-free  state  is  approached.  The  small,  almost  stress- 
free  elements  then  no  longer  fit  together  to  form  a  continuous  body.  As  the 
size  of  the  little  elements  approaches  zero,  the  truly  stress-free  state  is 
approached.  Also,  as  the  size  of  the  elements  diminishes,  the  aforementioned 
geometrical  mismatch  loses  its  interpretation.  It  is  replaced  by  the 
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assumption,  however,  that  in  the  limit  of  elements  with  vanishing  dimension  the 

configuration  of  the  body  is  the  result  of  an  incompatible  deformation  (in  the 

2 

sense  of  integrability  of  the  deformation  tensor)  of  a  continuous  body.  This 
intermediate  state  is  called  the  reference  configuration  because  it  is  the 
(time- varying)  state  of  vanishing  elastic  deformation. 

In  an  attempt  to  develop  the  desired  relations  with  a  minimum  of  confusion, 
notational  rules  to  be  followed  in  selecting  indices  for  the  various  tensor 
kernels  are  introduced.  Upper  case  Latin  indices  relate  to  the  initial  config¬ 
uration,  lower  case  Greek  to  the  reference  configuration,  and  lower  case  Latin 
to  the  current  configuration. 

Let  Xa  be  the  coordinate  in  the  reference  configuration  of  the  material 
particle  which  was  originally  at  Xq.  The  deformation  from  the  initial  to  the 
reference  state  cannot  be  written  as  a  differentiable  function  of  material 
coordinate,  as  in  (2),  because  of  the  assumption  of  incompatibility.  Indeed, 
it  is  precisely  the  lack  of  such  a  deformation  function  that  is  implied  by 

incompatibility.  The  plastic  deformation  of  a  line  element  dXQ  emanating 

K  ot  ct 

from  Xq  in  the  initial  configuration  into  line  element  dX  at  X  is 

defined  by  the  deformation  tensor  °bat,  through 

dXK  =  °b“  dXK  (7) 

K  O 

Again,  because  of  the  assumption  of  incompatibility,  °bav.  is  not  a  matrix  of 
partial  derivatives,  regardless  of  the  coordinate  system  being  used. 

Since  the  plastic  deformation  carries  the  particle  with  material  coordinate 

K  a  £ 

Xq  to  the  place  X  ,  and  the  total  deformation  (2)  carries  it  to  x  ,  the 


2 

The  foregoing  argument  is  not  a  mathematical  one,  but  a  physical  one  leading 
to  mathematical  assumptions. 
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elastic  deformation  carries  it  from  Xa  to  x1.  This  is  a  deformation  from 
an  incompatible  configuration  to  the  compatible  current  state  and,  therefore, 
the  tensor  describing  this  deformation  is  again  not  a  matrix  of  partial  deriva 
tives,  regardless  of  the  coordinate  system.  The  elastic  deformation  of  the 
line  element  dXa  at  X°  into  the  element  dx1  at  x1  is  defined  by  the 
tensor  b1^  through  the  non-integrable  relation 

dx1  =  b1  dXa  (8) 

a 

Substituting  dXa  from  (7)  into  (8)  there  results 


,  x  .1  o,  a  K 
dx  =  b  b  „  dX 
a  K  o 


Comparing  this  with  (3)  and  recalling  the  continuity  assumption,  we  get  a 
relation  between  deformation  tensors 


°b° 


(9) 


Thus ,  even  though  the  elastic  and  plastic  parts  of  the  deformation  are  defined 
by  tensors  which,  in  general,  are  not  gradients,  the  total  deformation  is 
described  by  a  tensor,  each  component  of  which  is  a  partial  derivative.  A 
simple  physical  example  of  this  phenomenon  where  incompatible  deformations 
combine  to  yield  a  compatible  one  is  easily  devised.  Consider  a  plate  of 
elastic-plastic  material,  out  of  which  we  cut  a  piece  of  material,  say  of 
square  shape.  Loads  are  applied  to,  and  then  removed  from,  the  cut  out  piece, 
resulting  in  homogeneous  plastic  deformation  and  no  residual  elastic  deforma¬ 
tion.  The  plate  now  resides  in  a  configuration  resulting  from  nonhomogeneous 
plastic  deformation  only.  This  state  does  not  satisfy  the  continuity  axiom, 
however,  because  we  had  to  introduce  the  cut  to  achieve  it.  To  restore  the 
plate  to  a  continuous  configuration  the  cut  out  piece  is  welded  back  into  the 
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vacant  square,  so  that  initially  adjacent  particles  are  again  adjacent, 

assuming  this  can  be  done  without  causing  further  plastic  deformation.  A 

residual  elastic  deformation  will,  in  general,  result.  With  this  example  in 

mind,  an  interpretation  of  (9)  is  that,  even  though  nonhomogeneous  plastic 

deformation  alone  leaves  a  body  of  elastic-plastic  material  in  a  physically 

impossible  configuration,  the  body  takes  advantage  of  its  elasticity  to  relieve 

this  state  so  that  it  can  reside  in  a  continuous  configuration. 

The  Green's  tensors,  which  are  the  metric  tensors  for  the  various  deformed 

configurations,  may  now  be  written  in  terms  of  the  deformation  tensors  a"^, 

b1  ,  °b®.  Consider  first  the  total  deformation.  To  express  the  current 
a  K 

length  of  a  line  element  ds  in  terms  of  the  original  element  in  the  initial 
configuration  dXQ,  (3)  is  substituted  into  (4),  yielding 

ds2  =  CKL  dXo  dXo  *  CKL  =  gija  Ka^L  *  (10) 

It  is  clear  from  (10)  that  is  a  metric  tensor  on  our  original  point  space. 

Furthermore,  it  can  be  shown  to  be  a  Euclidean  metric,  that  is,  to  have  a 

i  L  i  K 

vanishing  curvature  tensor,  by  using  the  fact  that  3a  K/3Xo  =  <*a  ^/3X  '.  In 
terms  of  material  coordinates,  the  total  strain  E^>  defined  by 

ds2  -  dS2  =  2E  dXK  dXL  ,  (11) 

o  KL  o  o 

may  be  obtained  from  (1)  and  (10)  as 


2EKL  =  CKL  '  oGKL 


(12) 


As  will  subsequently  be  seen,  the  total  strain  cannot,  without  special  inter¬ 
pretation,  be  written  as  the  sum  of  elastic  and  plastic  strain,  as  in  the  case 
of  infinitesimal  deformation.  In  a  similar  manner,  the  current  length  of  a  line 
element  in  terms  of  the  reference  configuration  is  given  by 
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ds2  =  B  0dXadxe  ,  B  Q  =  g^b1  b1.  ,  (13) 

ctg  ag  a  g 


and  the  length  of  a  line  element  in  the  reference  configuration  dS  in  terms 
of  the  initial  configuration  is 


dS2  =  B..TdXK  dXL  ,  Bv. 

o  KL  o  o  ’  o  KL 


o.o  o  g 

Gog  b  K  b  L  ’ 


K  L  K 

where  G  Q  =  G  g  g  0  and  g  is  the  appropriate  Euclidean  shifter.  The 
ap  o  KL  a  p  a 

two  metric  tensors  are,  in  general,  non-Euclidean  metrics.  The  corresponding 
elastic  and  plastic  strains  are  given,  respectively,  by 


Bag  ‘  Gag  * 


oBKL  "  oGKL 


The  fact  that  the  elastic  and  plastic  strains  are  not  additive  is  readily 
demonstrated.  By  using  relations  (9),  (10),  (13),  (14),  (15)  it  can  be  shown 


CKL  =  „»KL  *  VK  “A  ' 


Subtracting  QGKL  from  both  sides  and  making  the  substitution,  which 

jfe  (j 

defines  b 


o,a  _  a  ^  *  a 
b  K  ~  gK+  b  K 


we  obtain 


p  pe  a  g  e  ,a  *  g  g  *  a  *  a  *  g  . 
EKL  +  Eagg  Kg  L  +  Eag(g  KbL+gLbK+  b  K  b  L 


The  left  side  of  (16)  is  the  total  strain  in  terms  of  the  material  coordinates. 
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The  first  and  second  terms  on  the  right  side  are  the  plastic  and  elastic 
strains,  respectively,  in  terms  of  material  coordinates.  An  additional  term 
appears  on  the  right  side  because  initially  the  elastic  strain  is  measured 
with  respect  to  the  varying  reference  configuration.  Accordingly,  the  trans¬ 
formations  indicated  must  be  made  to  express  elastic  strains  in  terms  of 
material  coordinates.  In  the  case  of  infinitesimal  deformation,  this  additional 
term  is  of  "higher  order"  than  the  others,  the  distinction  between  various  con¬ 
figurations  becomes  negligible,  and  the  total  strain  is  simply  the  sum  of 
elastic  and  plastic  strains.  Green  and  Naghdi  [11]  write  the  finite  strain 
tensor  as  the  sum  of  an  elastic  and  a  plastic  part.  Their  elastic  strain  is 
different  from  ours,  that  is,  it  must  include  all  but  the  first  term  on  the 
right  side  of  (16)  and  therefore  must  depend  on  the  plastic  deformation  (as 
Green  and  Naghdi  subsequently  assume). 

Consideration  thus  far  has  been  limited  to  the  geometrical  properties  of 
the  deformation.  Attention  is  now  directed  toward  description  of  the  kinematic 
features.  The  fundamental  measure  of  rate  of  change  of  configuration  will  be 
taken  to  be  the  stretching  tensor,  which  defines  a  measure  of  the  rate  of 
change  of  current  length  per  unit  current  length.  This  tensor  then  defines  a 
measure  of  instantaneous  rates ,  and  thus  depends  only  on  the  current 
configuration. 

To  proceed  it  is  necessary  to  introduce  a  few  new  tensors,  related  to 

K 

those  previously  defined.  These  are  the  inverse  deformation  tensors  p 
a  o  K 

q  i ,  q  which  are  defined  by  the  systems  of  linear  equations 


L 

P 


K 


n  l  a 


o  L  o,  a 

q  b  v 

^  a  K 


(17) 


and  which  define  the  locally  affine  transformations 
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dX 


=  p  .dx 
r  1 


dX 


8  _ 


8  .  i 

q  .dx  , 


dX 


°qL  dX° 
^  a 


(18) 


These  inverse  deformation  tensors  are  obtained  by  solving  (17)  by  Cramer's 
rule;  for  example, 

p  .  =  cof(a1  )/det(a-1  )  .  (19) 

X  L  1\ 

Also,  observation  of  the  fact  that  the  metric  tensors  do  not  depend  explicitly 
on  time  and  have  vanishing  covariant  derivatives  allows  us  to  treat  the  metric 
tensors  as  constants  when  performing  material  differentiation.  The  vanishing 
of  the  covariant  derivative  of  the  metric  tensor  can  be  proven  by  direct  sub¬ 
stitution  or  by  observing  its  truth  in  the  underlying  Cartesian  coordinate 
system  and  invoking  the  invariance  property  of  tensor  equations. 

The  tensor  D..  ^  is  called  the  stretching  tensor  and  is  defined  by 
• 

ds2  =  2D^dx*dx^  ,  (20) 

where  the  dot  denotes  material  derivative.  Substitution  of  (10)  into  (20) 
and  performance  of  the  indicated  differentiation  leads  to 


11  K  L 
2  CKLP  iP  j 


t  -k  K  . 
sym(gika  Kp  . ) 


(21) 


The  stretching  is  thus  linear  in  the  material  derivative  of  the  deformation 
tensor  which  implies  that  it  is  linear  in  the  gradient  of  the  velocity  and, 
therefore,  linear  in  the  velocity  itself.  We  cannot,  however,  expect  to  obtain 
a  split  of  the  stretching  into  a  sum  of  elastic  and  plastic  parts  by  writing 
the  total  particle  velocity  as  a  sum  of  the  two  particle  velocities  (which 
do  indeed  exist)  corresponding  to  the  elastic  and  plastic  parts  of  the 
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deformation,  because  these  two  velocities  cannot  be  written  as  differentiable 
fields.  The  assumption,  tacit  in  the  statement  following  (21),  of  the  exis¬ 
tence  of  the  usual  kinematic  relation  stating  that  the  time  derivative  of  the 
deformation  tensor  equals  the  gradient  of  the  velocity  is  invalid  for  each  of 
the  elastic  and  plastic  deformations  separately. 

The  total  stretching  tensor  can,  however,  be  separated  into  the  sum  of 
an  elastic  stretching  tensor  and  a  plastic  stretching  tensor.  Substitution  of 
(9)  into  (21)  and  using  relations  between  metrics,  such  as  (6),  leads  to  the 
result 


,  /k  a  .  ^  ,  ,  k  o'a  o  L  g  .. 

sym(gikb  aq  +  sym(g.kb  a  b  L  g  ^q  . ) 


,  ,  k  a  ..  y  \>  a  ,k  o*g  o  L  y  m  .. 

sym(g..b  q  .)  +  g  .g  .  sym(G  g  ,b  .  b  ,  q  q  g  )•  (22) 

3  6ik  ]  6  Is  3  3  ave  kg  L  H  y  m6  p 


It  is  clear  that  the  first  term  on  the  right  side  of  (22)  defines  the  instan¬ 
taneous  rate  of  deformation  of  the  current  configuration  with  respect  to  the 
reference  configuration,  while  the  second  term  defines  the  instantaneous  rate 
of  deformation  of  the  current  configuration  due  to  the  changing  reference  con¬ 
figuration.  The  former  term  is  therefore  the  elastic  stretching  and  the  latter 
is  the  plastic  stretching.  The  plastic  stretching  takes  on  a  much  more  compli¬ 
cated  form  than  the  elastic  because  it  is  actually  the  part  of  the  instantaneous 
motion  of  the  current  configuration  in  terms  of  the  instantaneous  motion  of  the 
reference  configuration.  The  interior  term  °b^L°qL^  characterizes  the  motion 
of  the  reference  configuration.  This  term  is  then  contracted  with  b  .  and 

p 

qYm  which  deforms  the  characterization  into  the  current  configuration,  where 
the  actual  measurement  is  made.  Representing  the  elastic  and  plastic  stretch- 

0  D 

ings  by  D..  and  D. .,  (22)  is  written  as 
°  rj 


D.  .  =  D*.  +  D?. 

13  13  13 


D?.  +  g".gV.DP 
i&  3  PV 


(23) 
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One  additional  quantity  of  interest  is  the  stretching  of  the  reference 

s 

configuration  (i.e.,  the  measure  of  plastic  strain  rate).  This  rate  D  „  is 

ap 

defined  by 

dS2  =  2DS0dXadXP  ,  (24) 

af3 

and  is  easily  determined  to  be 

D®  =  sym(G  °bYT°qL0)  •  (25) 

a0  J  ay  L  ^  g 

Following  Truesdell  and  Toupin,  this  rate  is  called  the  slippage  tensor  and  is 
different  from  the  plastic  stretching  above.  (The  terminology  probably  stems 
from  the  fact  that  plastic  flow  is  the  macroscopic  manifestation  of  progressive 
microscopic  slip,  or  dislocation  motion,  on  discrete  planes  in  the  crystal 
lattice. ) 

In  summary,  expressions  for  the  following  rates  have  been  derived: 

D? .  -  measure  of  instantaneous  motion  of  the  current 

^  configuration  due  to  elastic  deformation 

D?.  -  measure  of  instantaneous  motion  of  the  current 

^  configuration  due  to  plastic  deformation 

-  measure  of  instantaneous  motion  of  the  current 
configuration  due  to  plastic  deformation,  referred 
to  the  intermediate  configuration 

D  -  measure  of  instantaneous  motion  of  the  reference 
configuration . 

Furthermore,  although  only  stretching  of  lineal  elements  has  been  considered 
above,  shearing  (i.e.,  instantaneous  shear  rate)  can  also  be  determined  from 
these  tensors.  For  example,  let  ip  be  the  angle  between  elements  dx^  and 
dx*  in  the  coordinate  directions  at  x1.  For  any  two  elements 

ds,dsr.  cos  tl;  =  g..dx^  dx^. 

12  T  1  2 
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Then  at  any  instant  tji  =  ir/2  and 


♦ 


where  the  element  lengths  are  measured  as  in  (4). 


(26) 


RESTRICTIONS  IMPOSED  BY  PHYSICAL  ASSUMPTIONS 

In  the  preceding  section  the  geometric  and  kinematic  features  of  the 
deformation  of  a  body  of  elastic-plastic  material,  undergoing  finite  displace¬ 
ments,  were  considered,  without  regard  for  actual  deformation  processes  of 
such  material.  We  now  consider  some  of  the  simplifications  introduced  by 
applying  assumptions  resulting  from  experimental  observation. 

The  most  frequent  assumption  made  concerning  elastic-plastic  deformation 
is  that  the  volume  change  per  unit  volume  is  recoverable,  that  is,  the  volume 
change  due  to  plastic  deformation  is  zero.  To  determine  the  mathematical  form 
of  this  assumption,  consider  a  small  volume  element  of  material  dV  in  the 
reference  configuration  which  occupied  the  element  dVQ  in  the  initial  con¬ 
figuration.  In  general,  these  volume  elements  are  related  by 


dV  =  J  dV 
o  o 


(27a) 


where 


J 

o 


det  (G  )  1/2 

<°bV 

O  KJj 


(27b) 


To  satisfy  the  condition  of  no  plastic  volume  change,  set  J  =  1.  It  can 
also  be  shown  that 


J  2 
o 


KT 

III  =  det  (  G  BTlI) 
o  o  o  LM 
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so  that  an  equivalent  condition  is  IIIo  =  1*  Yet  another  equivalent  condition 

is  tr(GaYDSQ)  =  0. 

YB 

A  second  assumption,  concerning  the  elastic  or  recoverable  part  of  the 
deformation,  might  be  made.  As  stated  above,  the  material  volume  change  during 
elastic-plastic  deformation  is  completely  elastic,  and  we  can  conceive  of 
deformation  processes  during  which  this  dilatancy  is  large.  This  is  not  the 
case,  however,  for  elastic  shears.  During  elastic-plastic  deformation  of  most 
metals  yielding  begins  long  before  the  elastic  shears  increase  to  the  point 
where  they  must  be  considered  large.  It  seems  that  in  many  applications  the 
assumption  that  the  elastic  shears  remain  small,  even  though  the  plastic 
shears  become  large,  is  a  reasonable  one. 

Before  proceeding,  some  sort  of  meaning  should  be  attached  to  the 
description  "small."  Here  we  say  that  a  deformation  tensor  h1  is  small  if 


h 


max  ,i 
.  n 
i,a  a 


and  h  «  1  . 


The  above  assumption  is  given  a  mathematical  form  by  stating  that  the 
elastic  deformation  tensor  is  a  small  deviation  from  a  uniform  dilation.  This 
is  written  in  terms  of  physical  components  of  tensors  (denoted  by  hat  *  ) 
because  Bfi"^  are  the  physical  components  of  a  deformation  tensor  for  uniform 
dilation,  where  B  is  a  scalar  function,  0  <  B  <  00 .  The  assumption  takes 
the  form 

61  =  B(61  +  h1  )  ,  (28) 

a  a  a 

i  -i  3 

where  h  is  small.  Without  loss  of  generality,  we  can  take  tr(h  q)  =  0. 

3  •  •  i 

Note  that  the  corresponding  quantity  tr(h  )  cannot  be  formed. 

^Physical  components  of  a  deformation  tensor  measure  change  in  length  per  unit 
length,  while  tensor  components  measure  change  in  generalized  coordinate  per 
unit  generalized  coordinate  [10], 


16 


For  assume  this  is  not  true.  Setting 


B^CB  +  |  trU^lTh^  +  i  6^  tHh^)  , 


*  1  -i 

3  =  B  +  itrth^)  , 


relation  (28)  reduces  to 


A1  &  1  &  *  1 

b1  =  8  (61  +  h1  ) 

a  a  a 


where  h  is  small.  Redefining  the  starred  quantities  to  be  their  unstarred 

counterparts,  we  again  arrive  at  the  form  (28),  with  the  trace  of  h1^  vanish- 

2. 

ing.  The  nine  components  of  b  are  thus  replaced  by  8  and  the  eight 

*  i 

independent  components  of  h 

Reverting  to  tensor  components,  (28)  becomes 


bl„  =  6<thV  • 


To  determine  the  volume  change  of  the  deformation,  suppose  that  the  material 
volume  element  dv  in  the  current  configuration  occupied  the  element  dVQ  in 
the  initial  configuration;  then 


det(g..)  , 

dV  *  JdVo  •  J  ■  (dst (  "£■)>  det(aV  • 

O  KL 


Substituting  (9)  into  (30)  and  assuming  zero  plastic  volume  change,  J 
reduces  to 


det(g..)  , 

J  =  det(by>  • 

ap 
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Furthermore,  substituting  (29)  into  this  relation  and  expanding  in  a  nine¬ 
dimensional  Taylor  series  about  h  =  0,  we  obtain 

•  •  3  .  . 

det(b1  )  =  63{det(g1  )  +  £  h1  cofCg1  )}  +  o(h)  .  (31) 

l  ,a=l 

A  similar  expansion  about  h  =  0  yields 


cof (b1^)  =  32{cof(g1a)  +  ei;.keaeY  h3’ghkY}  +  o(h) 


(32) 


This  formula  is  of  use  in  solving  (17)  for 


by  Cramer's  rule. 


SPECIAL  DEFORMATIONS 

We  now  list  some  results  for  specific  coordinate  systems  and  special 

deformations.  In  particular,  plane  strain  in  Cartesian  coordinates  and  axially 

symmetric  strain  in  circular  cylindrical  coordinates  are  considered. 

In  Cartesian  coordinates  the  material  coordinates  of  an  arbitrary  particle 
K  12  3 

of  the  body  are  Xq  =  (Xq,  Xq,  Xq)  before  deformation  begins.  At  time  t 
this  particle  is  at  place  x1,  and  the  fact  that  the  deformation  is  plane 
strain  (assumed  independent  of  the  3  coordinate)  yields  the  functional 
relations 


1  l,vl  2  .  2  2.1  2 

x  =  x  (X  ,  X  ,  t)  ,  x  =  x  (X  ,  X  ,  t) 

O  o 7  0*0* 


The  deformation  tensor  has  the  form 


[a^]  = 


0 

0 


(33) 


0 


0 


1 
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The  elastic  and  plastic  deformation  tensors  b1  and  °ba  have  the  same 

Ot  K 

form  as  (33),  except  for  the  33  components.  These  are  not  each  unity  but  are 
related  by  b^3°b^3  =  1. 

The  metric  tensor  and  Euclidean  shifters  reduce  to  unit  tensors  for 
Cartesian  coordinate  systems. 

The  condition  that  the  plastic  volume  change  is  zero  can  be  stated, 
from  (27b),  for  this  case  as 


o,  3  ,o,l  o,2  o,l  o,  2  . 
b3(bl  b  2  ~  b  2  bl)  " 


Relation  (35)  may  be  used  to  eliminate  °b33  from  all  other  relations.  The 
condition  that  elastic  shears  are  small  becomes 


b1  =  0(61  +  h1  ) 

a  a  a 


(36a) 


where 


(36b) 


When  the  motion  is  prescribed  in  terms  of  material  coordinates  in  a 
Cartesian  system  the  material  time  derivative  reduces  to  the  usual  partial 
derivative  with  respect  to  time,  and  deformation  rates  may  be  easily  calculated. 
Also,  in  Cartesian  coordinates,  the  covariant  derivative  reduces  to  the  usual 
partial  derivative.  For  the  case  of  one-dimensional  strain  the  geometric  and 
kinematic  relations  are  given  by  Lee  and  Liu  [7], 

As  a  second  example  consider  a  deformation  which  is  rotationally  symmetric 
about  an  axis.  The  circular  cylindrical  coordinates  of  an  arbitrary  particle 
in  the  initial  configuration  are  =  (Rq,  0o,  2q),  and  this  particle  occupies 
the  place  x1  =  (r,  0,  z)  at  time  t,  with  the  Zo~axis  being  the  axis  of 
rotational  symmetry.  The  deformation  has  the  form 
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r  =  r(R  ,  Z  ,  t)  ,  6  =  0  ,  z  =  z(R  ,  Z  ,  t)  (37) 

o  o  o  o  o 

and 

[a1,]  = 

As  before,  b1  and  °ba  have  the  same  form  as  (38)  except  that  the  22 

Ot  I\ 

2  o  2 

components  are  not  each  unity,  but  that  b  2  °b  ^  =  1- 

The  metric  tensor  and  the  Euclidean  shifters  depend  on  position  in  this 
coordinate  system,  and  may  be  determined  as 


(38) 


CoGKL] 


(39a) 


where  i  and  K  are  not  summed. 

The  condition  that  there  be  no  volume  change  due  to  plastic  deformation  is 
determined  from  (27b)  to  be 

o,  2  ,o,l  o,3  o,l  o,  3  .  _  ,, 

b  2(  b  1  b  3  ~  b  3  b  1}  =  VR  (40) 

where  R  is  the  radial  coordinate  in  the  intermediate  configuration 
X°  =  (R,  0,  Z). 

In  this  coordinate  system,  the  covariant  derivative  is  written  in  terms  of 
the  Christoffel  symbols  .  Of  the  27  components  only  three  are  nonzero,  and 


these  are 
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CONCLUSION 

We  have  considered  the  finite  deformation  of  a  body  of  elastic-plastic 
material  by  introducing  a  particular  intermediate  configuration,  distinct  from 
the  initial  and  current  configurations  of  the  deforming  body,  which  is  the 
state  of  the  body  due  to  its  plastic  deformation  alone.  The  main  result  is 
that  the  total  rate  of  deformation  tensor,  in  an  orthogonal  curvilinear 
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coordinate  system,  can  be  written  as  the  sum  of  rates  of  elastic  and  plastic 
deformation  as  shown  in  (22).  The  elastic  stratching  takes  on  an  expected 
form,  while  the  plastic  stretching  tensor  is  not  the  intuitive  rate  of  defor¬ 
mation  (which  is  called  the  slippage  tensor  above).  That  is,  the  plastic 
stretching  is  not  merely  a  shift  to  the  current  configuration  of  the  rate  of 
deformation  of  the  reference  configuration.  It  might  be  called  instead  a 
deformation  of  this  rate,  which  reduces  to  a  simple  shift  only  when  b1^  =  g1Q» 
Whereas  a  shift  leaves  the  proper  numbers  and  proper  vectors  of  a  tensor  un¬ 
altered,  the  proper  numbers  and  vectors  of  are  different  from  those  of 

s  »  ,  s  a  B 

D  0  and,  in  general,  D%.  t  D  0g  .g  .. 

aB  ’  6  i]  aB  i  D 

Having  established  the  additivity  of  the  rates  of  elastic  and  plastic 
deformation,  the  general  scheme  for  development  of  proper  forms  for  constitu¬ 
tive  equations  becomes  clear.  The  particular  rates  of  (22)  must  be  written  in 
terms  of  material  properties,  stress,  and  temperature,  and  then  added  to  yield 
a  stress-deformation-temperature  relation  for  the  material.  In  a  recent  paper, 
Perzyna  [12]  assumed  an  additivity  resembling  that  of  (22)  and  from  there  pro¬ 
ceeded  to  discuss  forms  for  constitutive  equations  for  viscoplastic  materials. 
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